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Abstract 



We propose a supersymmetric A4 x SU(5) model of quasi-degenerate neutrinos which predicts 
the effective neutrino mass m ee relevant for neutrinoless double beta decay to be proportional 
to the neutrino mass scale, thereby allowing its determination approximately independently 
w—> of unknown Majorana phases. Such a natural quasi-degeneracy is achieved by using A4 family 

^~H, symmetry (as an example of a non-Abelian family symmetry with real triplet representations) to 

enforce a contribution to the neutrino mass matrix proportional to the identity. Tri-bimaximal 
neutrino mixing as well as quark CP violation with a « 90° and a leptonic CP phase <5mns ~ 90° 
arise from the breaking of the A4 family symmetry by the vacuum expectation values of four 
"Havon" fields pointing in specific postulated directions in flavour space. 



E-mail: antusch@mppmu.mpg.de 
E-mail: king@soton.ac.uk 
E-mail: spinrath@mppmu.mpg.de 



1 Introduction 

Over the past dozen years or so our knowledge of the neutrino sector has increased dramatically 
with the discovery of neutrino mass and mixing in atmospheric and solar neutrino oscillations, 
followed by the observation of terrestrial neutrino oscillations in long baseline neutrino experiments 
which have confirmed and refined the earlier results []J. Yet, despite this progress, which may even 
be termed a "neutrino revolution", there are many questions about neutrinos which remain unan- 
swered. Perhaps the most pressing of these is the origin, nature and magnitude of neutrino mass, 
since neutrino oscillations only provide information about the squared mass differences between 
neutrino species which are independent of the absolute neutrino mass scale or the nature of the 
neutrino mass (i.e. Dirac or Majorana). In the absence of any confirmed experimental signal from 
either beta decay end-point experiments or neutrinoless double beta decay experiments, the most 
stringent limits on the absolute neutrino mass scale comes indirectly from cosmology where one 
typically obtains the limit on the absolute neutrino mass scale expressed in terms of the lightest 
neutrino mass as m^ghtest ^S 0.2 eV [2J. Thus, there remains the interesting possibility that neutri- 
nos are quasi-degenerate, which one may roughly define as miightest > 0.05 eV, where the lower limit 
is approximately set equal to the square root of the atmospheric neutrino mass squared difference. 

The current generation of running or planned neutrinoless double beta decay experiments are 
capable of discovering quasi-degenerate neutrinos, as defined above, within the next years. Such 
a discovery would herald a new neutrino revolution to rival the last one, and would lead to an 
explosion of interest in theoretical models capable of accounting for quasi-degenerate neutrinos. In 
general having quasi-degenerate neutrinos does not lead to a sharp prediction for the neutrinoless 
double beta decay observable m ee as a function of mi ~ m<i ~ m.3 due to the presence of unknown 
phases in the neutrino mass matrix [3|]. The general conclusion that unknown phases enter the 
prediction for m ee also remains valid in models which combine the experimental observation of (at 
least approximate) tri-bimaximal (TB) [J] lepton mixing with the possibility of a quasi-degenerate 
neutrino mass spectrum. The reason is that the relevant phases entering m ee are Majorana phases. 
Allowing for arbitrary Majorana phases and considering a quasi-degenerate neutrino mass spectrum 
and TB mixing, m ee can still be in the approximate interval m ee 6 [^lightest /3, ^lightest]- 

What would we learn about the origin of neutrino mass from the discovery of quasi-degenerate 
neutrinos in neutrinoless double beta decay? Clearly this would imply that neutrinos are Majorana, 
and possibly (but not necessarily) that would indicate that a seesaw mechanism is at work, but 
what kind of seesaw mechanism, i.e. is it type I or II'^J There are known examples of type I and 
type II seesaw models which can lead to quasi-degenerate neutrinos as well as TB lepton mixing, 
so clearly quasi-degenerate neutrinos would not distinguish different types of seesaw mechanism. 
For example, the super symmetric (SUSY) Grand Unified Theory (GUT) based on 50(10) with 
family symmetry PSL(2, 7) proposed in [5J is based on the type II seesaw mechanism and leads 
to TB mixing and allows quasi-degenerate neutrinos. On the other hand the SUSY A4 model in 
[6|] based on the type I seesaw mechanism also leads to TB mixing and allows quasi-degenerate 
neutrinos. Interestingly, the SUSY A± x SU(5) model with type I seesaw mechanism does not 
favour quasi-degenerate neutrinos [7j], whereas a related model with a type II seesaw mechanism 
does allow quasi-degenerate neutrinos [8j|. More generally, there is a huge literature on family 
symmetry models based on A4 [9(] or other symmetries [id ]. However, to our knowledge, in all 
above examples, quasi-degenerate neutrinos are subject to the mentioned phase uncertainties in 
the prediction of m ee as a function of mi. 



1 We shall not consider the type III or further types of seesaw mechanism in this paper. 



It is interesting to ask, in what class of theories would we learn the most about the neutrino 
mass scale miightest from the discovery of a measurement of m ee in neutrinoless double beta decay? 
Clearly the answer would be those theories which predict m ee uniquely as a function of mi ig htest 
without ambiguities from unknown phases, but the next question is do such theories exist? Perhaps 
surprisingly the answer is in the affirmative, and, even more surprisingly, the class of theories which 
have this property turn out to suggest the way that the neutrino mass matrix is generated, namely 
by a usual type I seesaw contribution with two or three right-handed neutrinos, plus an additional 



contribution proportional to the unit matrix. In ll|, two of us proposed a class of theories of 
exactly this kind which we referred to as a "type II upgrade of type I seesaw models" . In this class 
of models the additional contribution to the neutrino mass matrix was realised by an additional 
type II seesaw. The type I seesaw part of the neutrino mass matrix, which controls the mass squared 



differences and mixing angles, was governed by sequential right-handed neutrino dominance 121 ] . 
The effect of such an additional unit matrix structure implies that for quasi-degenerate neutrino 
masses the Majorana CP phases are small and thus m ee ~ Slightest- Although the class of models 
was specified, no realistic "type II upgrade" model has ever been proposed. 

In this paper we shall propose a model following the idea of an additional contribution the 
neutrino mass matrix proportional to the unit matrix based on A4 family symmetry with SU(5) 
grand unification. The model contains tri-bimaximal neutrino mixing after the A4 family symmetry 
is broken as an indirect result of the assumed aligned "neutrino flavons" in the type I seesaw sector 



via constrained sequential dominance 13l |. These "neutrino flavons" break the A^ symmetry, being 



assumed to be aligned along the columns of the TB mixing matrix, but quadratic combinations 



of the neutrino flavons respect accidentally the neutrino flavour symmetry as discussed in [141 ]. 
Further "quark flavons" are assumed to be misaligned compared to the "neutrino flavons" and are, 
together with the "neutrino flavons", responsible for quark and charged lepton masses and quark 
mixings. As expected, due to a possibly large type II seesaw contribution, or alternatively due to 
an additional type I seesaw contribution from an additional triplet representation of right-handed 
neutrinos, the model can naturally predict the neutrinoless double beta decay mass observable to 
be approximately equal to the neutrino mass scale. We also make a detailed fit to quark masses 
and mixing using the misaligned quark flavons and show that a simple ansatz for the phase of 
one of the misaligned quark flavons leads to successful quark CP violation. In order for radiative 
corrections not to modify too much the TB mixing for quasi-degenerate neutrinos [l5l |. we shall 
restrict ourselves to low values of the ratio of Higgs vacuum expectation values (vevs) tan/3 < 1.5. 
For such low tan/3 < 1.5, a viable GUT scale ratio of y^/y s is achieved within SUSY SU(5) GUTs 
using a Clebsch factor of 9/2, as proposed recently by two of us in [16(. For the third generation 
we use b-r Yukawa coupling unification y T /yb = 1 at the GUT scale which is viable for low tan/3 



(see, e.g., [17]]). 

The layout of the remainder of the paper is as follows. In section 2 we present the model. In 
section 3 we perform a numerical fit to the quark and charged lepton masses and quark mixing 
angles and CP violating phase and discuss the neutrino masses and lepton mixing angles. Section 
4 summarises and concludes the paper. In Appendix A we give a renormalisable superpotential 
and explicit expressions for the effective couplings and Appendix B contains a possible vacuum 
alignment. 
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Table 1: Representations and charges of the superfields. The subscript i on the fields Ti, iV; and d is a family index. 



The flavon fields 



i>23 can be associated to a family via their charges under . 



I4. The subscripts on the Higgs 



fields H, H and extra vector- like matter fields A, A denote the transformation properties under SU(5). 

2 The Model 

In this section, we propose and describe a SUSY GUT model based on the unified SU(5) gauge 
group as well as on the family symmetry A4 amended by some discrete Z| x Z| symmetries and an 
U(1)r symmetry as specified in Tab. [TJ 



2.1 Symmetries and Field Content of the SU(5) GUT Model 

Let us start introducing the model by specifying the field content and the symmetries. The Standard 
Model matter fields fit nicely into the two representations 5, which we call F, and 10, which we 
call T. Explicitly they are given as 
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(2.1) 



where the lower indices R, B and G denote the quark colours and i = 1,2,3 is the family index. 
In our model, we consider that the three generations Fi form a triplet representation 3 of an A4 
family symmetry whereas the three generations T, form singlets 1 under A4G In the following, we 
suppress the A4 indices. In addition, we consider two right-handed neutrinos, singlets under SU(5) 
as well as under A4, labeled by N\ and N2. 

We furthermore consider fifteen-dimensional Higgs representations -Hi5, H15 which contain 
SU (2) L-tiiplet Higgs fields that obtain induced vevs after electroweak symmetry breaking. H\§ 
induces in this way a type II seesaw contribution to the neutrino mass matrix which is, to leading 
order, proportional to the unit matrix and can increase the neutrino mass scale without modifying 
the values for the leptonic mixing angles. 

SU(5) is spontaneously broken by the vev of the 4>23 field, electroweak symmetry is broken by 
the vevs of the Higgs fields H5, H5, H45, H^ and A4 is spontaneously broken by the vevs of the 
flavon fields, i.e. the family symmetry breaking Higgs fields </>i23, 023, 03 and 023- We comment 
below on the specific directions in which we assume A4 to be broken by the flavons. 

On top of that, we consider additional "messenger" fields which are heavy and which, after 
effectively integrating them out of the theory, give rise to higher-dimensional operators generating 
the Yukawa coupling matrices as well as the mass matrix of the gauge singlet (right-handed) 
neutrinos iVj. 

The field content of our model as well as the symmetries are specified in Tab. [TJ We note 
that it is always possible to replace any product of commuting discrete symmetries by a single 
Abelian group U(l) with a suitable choice of charges for the fields, so it is possible to replace the 
Z| x Z| symmetry by a single U(l) symmetry, with an appropriate choice of charges. Indeed many 
models in the literature use an Abelian U(l) symmetry rather than a product of Zjy symmetries to 
control the operators. Although this looks simpler, it should be remarked that firstly an Abelian 
symmetry has infinitely many more group elements than any discrete symmetry, and secondly one 
must then confront the question of Goldstone bosons once the assumed global Abelian symmetry 
is broken. If the Abelian symmetry is gauged one must further complicate the model by ensuring 
that it is anomaly free. Therefore an auxiliary discrete symmetry, even a large one, has definite 
advantages over an Abelian symmetry. Furthermore discrete symmetries are ubiquitous in string 
theory constructions. Finally, the auxiliary discrete symmetry used here is rather a simple one 
consisting of a product of Z2 and Z4 parity factors. Thus we regard the use of the discrete Z2 and 
Z4 symmetries as being a well motivated, simple and attractive alternative to the use of an Abelian 
U(l) symmetry. 

We would like to remark that we do not explicitly consider the full flavour and GUT Higgs 
sector of the model and just assume that the SU(5) and A<± breaking vevs are aligned in the desired 
directions of field space. We assume that in these sectors issues like doublet-triplet splitting are 
resolved. Without specifying these sectors, a reliable calculation of the proton decay rate must also 
be beyond the scope of the present paper. The focus of the present paper is thus to illustrate that 
quasi-degenerate light neutrino masses can be realised together with a type II seesaw in a SU(5) 
GUT framework. 

We would furthermore like to remark that in addition to the type II seesaw contribution there 
is a possible additional contribution to the neutrino mass matrix proportional to the unit matrix 
from the messenger field A\ which is a singlet under SU(5) and a triplet under A\. When it is 



2 We note that in principle any non-Abelian family symmetry with real triplet representations, like, e.g., SO(3), 
would in principle be suitable for the construction of models with additional contributions to the neutrino mass 
matrix proportional to the unit matrix. In this paper we focus on A4, as a specific example. 
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Figure 1: Supergraph diagrams inducing the effective superpotential operators for the down- type quarks and charged 
leptons. 

integrated out, it also induces a contribution to the neutrino mass operator which is proportional 
to the unit matrix. 

2.2 The Effective A 4 x SU(5) Symmetric Superpotential 

The renormalisable superpotential resulting from Tab. [I] is given in the Appendix. Integrating out 
the heavy messenger superfields denoted by A, B and C, the Feynman diagrams in Figs. [HE] and [3] 
then lead to the effective non-renormalisable superpotential terms in the SU(5) and A4 unbroken 
phase: 
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After GUT symmetry breaking the SU(2)i doublet components from H§ and I/45 respectively H§ 
and i/45 mix and only the light states acquire the SU(2)i breaking vevs which give the fermion 
masses, as discussed in Appendix A where the effective couplings a appearing in the effective 
superpotential are also explicitly given. 
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Figure 2: Supergraph diagrams inducing the effective superpotential operators for the up-type quarks. 

2.3 Assumed Vacuum Alignment 

In the following we assume that the vevs of the A4 breaking flavon fields point in the following 
directions in field space such that 
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The bi are defined in Appendix A. These relations also define the quantities £123, £23 and €3. The 
breaking of A4 along the field directions of ^123 and (^93 allows us to realise tri-bimaximal neutrino 
mixing via constrained sequential dominance (CSD) [12]. It is also worth noting that the flavon 
vevs (^123) and (^23) are orthogonal, causing some of the terms in the superpotential to give a 
vanishing contribution to the mass matrices. In the following, we assume that CP is only broken 
spontaneously by the vev of the flavon ^23 • 

For the flavon ^23 one may suppose a priori a less constrained alignment, 
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(2.9) 



However empirically we find that the numerical fit to quark masses and mixings, in particular quark 
CP violation, seems strongly to prefer that the vacuum alignment of the flavon (^23 has its second 
component along the imaginary direction. To simplify the results of the numerical fit we shall 
restrict ourselves to the case: 

v = -i . (2.10) 
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Figure 3: Supergraph diagrams inducing the effective superpotential operators for the neutrino sector. 



In some future more ambitious theory one may attempt to reproduce Eq. (|2.10p as a result of some 
special vacuum alignment, but here we shall simply regard it as a special choice, or ansatz, which 
leads to a successful fit to quark CP violation. 

In Appendix B, we will discuss another possibility, namely to realise the flavon ^23 effectively 
by splitting it up into two flavons fa and fa, where one gets a purely real and the other a purely 
imaginary vev. For the effective superpotential in Eqs. (|2.2|) - (|2.7[) this would correspond to simply 
replacing ^3 — > <\>\ + 4>\ an d ^23-^45 — > faH^ + faH' 5 , with an additional Higgs field H' 5 . The field 
content of this extended version of the model that now includes also a vacuum alignment sector, is 
presented in Appendix B in Tab. 01 The predictions of the two model variants are identical at the 
level of precision discussed here. 

3 Numerical Fit to Fermion Masses and Mixings 

3.1 The Quark and Charged Lepton Sector 

We define our conventions for the Yukawa matrices such that the operators of the form FT(j)H and 
T 2 (j) 2 H give the following Yukawa terms in the Lagrangian: 

£y u k = -(. Y d)ijQidjH d - (J*)ijLiejH d - {Y*)ijQiUjR u + h.c. , (3.1) 

where the SU(5) relation Y d = Y^F would be fulfilled, if all Clebsch-Gordan factors were one. The 
convention we use here is the same as the one used by the Particle Data Group [33]. 

From Eqs. $T2$), ([23]), (pl8jh ([2"U]) and (f2T0j) the Yukawa matrix coupling the up-type quarks 
to the light up-type Higgs doublet with the b coefficients as defined in Appendix A is given as 

/2&ue| 3 h 3 e 23 e 3 \ 

Y u =\ 3& 22 e? 23 + (w 2 - l)b 2 2 e 2 3 b 23 e 123 e 3 , (3.2) 

\fri3e23e3 &23 e 123£3 ^33 / 

whereas the Yukawa matrices coupling the down-type quarks and charged leptons to the light 
down-type Higgs doublet are given as 

/ e 23 -£23 \ 

Y d = £123 £123 + 1623 £123 + W€ 23 , (3.3) 

\ e 3 / 

(0 C 2 3£23 -C23E23 \ 

C123C123 Ci23£l23 +iC23£23 Ci23£l23 + ^C 2 3£23 , (3.4) 

C 3 £ 3 / 

where c 3 , on, c 23 and C123 are the Clebsch-Gordan factors arising from GUT symmetry breaking, 
see, e.g., [la]. We have used the orthogonality of (^23) an d (^123) and considered the above 
described notation for the flavon vevs. We note that in the definition for the Yukawa matrices we 
have introduced a complex conjugation why here appears a phase factor of +i in the 2-2 elements 
of the down- type quark and charged lepton Yukawa matrices. 

With the given representations of the flavons, we obtain the following Clebsch-Gordan coeffi- 
cients 

C123 = 1 , c 23 = 1 , c 3 = 1 , C23 = 9/2 . (3.5) 
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Table 2: The model parameters for tan/? = 1.4 and Max 



500 GeV from a fit to the experimental data. 



For small values of tan/3 as we consider, the 1-loop SUSY threshold corrections are small and, 
taking the actual experimental values of the strange quark and muon masses into account, the 
GUT scale value of y^,/y s prefers C23 = 9/2, as argued in [16|] (see also [2QJ). 

From the charged lepton Yukawa matrix we can derive the following approximate relations for 
the eigenvalues 



lh 



C3£3 , V/j, = |C123£123 + i^^l 



lie 



C23€23Cl23£l23 



(3.6) 



Furthermore, since there is no 1-2 mixing from the up-sector, the mixing angle #12 is approximately 
given as 
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From those four equations the four e's can be calculated and the relation for the CKM phase gives 
at the GUT scale 

£23 



tan(5 C KM| 



£123 



1.22 



(3.8) 



The RG evolution of the measured value for 5ckm gives a GUT scale value of 1.20. So the value for 
the CKM phase (based on our assumed vacuum alignment) is already remarkably good if we only 
take the lepton masses and the value for #12 into account which are measured to high accuracyo 

For the detailed fit of the model to the data we applied the following procedure: We have taken 
the GUT scale Yukawa matrices from Eqs. (|3.2p . f|3.3[) and (|3.4p and calculated their RG evolution 
down to the scale mt(mt) for tan/3 = 1.4j and Msugy = 500 GeV with the REAP software package 
[261 ] . At the low scale we performed a x 2 fit to the quark masses and mixing and charged lepton 



We would like to remark that with the assumed spontaneous CP violation, real det Y u and det Yd and with the 
small I £23 1 = O(10 -4 ), the model might also provide a solution to the strong CP problem, along the lines discussed 



in [2 

4 



We note that in the MSSM small values of tan/3 are somewhat constrained due to bounds on the Higgs mass. 
However, we emphasize that our model may well be formulated in the context of the NMSSM or other non-minimal 
SUSY models where tan/3 of order one can readily be realised without these constraints. 
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Table 3: Fit results for the quark Yukawa couplings and mixing and the charged lepton Yukawa couplings at low 
energy compared to experimental data. A pictorial representation of the agreement between our fit and experiment 
can be also found in Fig. [4] 



masses depending on the parameters of the GUT scale Yukawa matrices. The fit gave a total y 2 
of about 3.5 where we have assumed a relative error of 1 % for the charged lepton masses and for 
the other observables we have taken the experimental errors. Since we have 11 parameters and 13 
observables this corresponds to a x 2 /dof of about 1.6. This is a good fit since we have neglected 
theoretical uncertainties like, e.g., threshold corrections which could be treated as additional errors 
on the data lowering the total x 2 . 

The results for the GUT scale parameters are listed in Tab.[2j We would like to remark that these 
parameters depend on tan/3 and Mstjsy and also are subject to several theoretical uncertainties. 
For example, we note that the Higgs fields H\§ and H\§ containing the Higgs triplets of the type 
II seesaw mechanism have masses at intermediate energy scale between Mqut and Mew- Their 
effects are not included in the RG analysis. The effects are small and may be neglected if j/a is 
small, but they could be sizable if y& is largeO Due to the additional theoretical uncertainties 
we do not explicitly give the errors on the high energy parameters or low energy fit results. The 
important input parameters for us are the charged lepton masses and quark mixing angles which 
have a experimental error much smaller than these uncertainties. 

In Tab. [3] the low energy results are shown and compared to experimental data. A graphical 
illustration is given in Fig. HI They illustrate that our minimal example model, with the assumed 
vacuum alignment of Eq. ()2.8p and (j2.9f) . can fit well the data. We turn now to the results for the 
neutrino sector. 



5 Since the coupling y/\ gives a contribution proportional to the unit matrix, it affects only the RG evolution of 
the mass eigenvalues, but not of the mixing angles. Nevertheless, the possibility of additional RG effects from j/a 
provides a theoretical uncertainty in our setup. 
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Figure 4: Pictorial representation of the deviation of our fit from low energy experimental data for the charged lepton 
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3.2 The Neutrino Sector 

The neutrino Yukawa matrix is obtained from Eq. (j2.4|) as 



^ 2 ei23\ 

Y v = ( b Ul e 2 3 b U2 ei2z ■ (3.9) 

-b Vl (-2Z b U2 ei23/ 

Additionally we have a diagonal mass matrix for the two right-handed neutrinos from Eq. (|2.7|) . 

and contribution proportional to the unit matrix coming from Eqs. (|2.5|) and (|2.6|) . 

(mo \ 
M x = m . (3.11) 

\ mo/ 

Using the seesaw relation 

m, = Mi - v^M^Yj , (3.12) 

we obtain for the neutrino mass matrix 

( l ° °\ m' Z 1 X l \ m 1 f° ° ( 
m„ = m 1 + -?- 1 1 1 + -^ 1 -1 | (3.13) 

\0 1/ \1 1 1/ \0 -1 1 

m = c v u + c v u —— , m 2 = -v u — and m 3 = -v v — . (3.14) 

^15 M4i b R2 b Rl 

In our model we therefore identify the neutrino masses as mi = mo, vti2 = mo + m^, m% = mo + m^, 
where without loss of generality we can take mo to be positive and real while m\, m% are real but can 
take either sign. With |mo| 3> \m{\, Im^l a quasi-degenerate mass spectrum of the light neutrinos 
can be explained in a natural way. In the following, we will mainly restrict ourselves to this case. 
From the structure of these matrices we see that we obtain tri-bimaximal mixing in the neutrino 
sector 



with 



^ 3 = , 0£ 3 = 45° , 9x2 = arcsin -= « 35.3° . (3.15) 



1 
73 



From the lepton sector we get the additional mixing contributions 

£123^123 



? 13 - ° ) #23 ~ ° i Fl2l 



C123C123 — 1C23 £ 23 



4.6° . (3.16) 



There is also a complex phase introduced by the charged lepton Yukawa matrix which can be 
calculated in the same way as in the quark sector 

5 e 12 = arctan C23£23 « -85.4° . (3.17) 

C123C123 
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27(: 



For the approximate calculation of the MNS mixing parameters at the GUT scale we can use 



MNS _ v ne 

s 23 ~ s 23 ~~ ^23 

s MNS e -< 3 NS ^ % _ ^e^ (3>18) 

where we have already discarded all trivial phases and RG corrections which we will discuss later. 
For the total leptonic mixing angles we obtain 
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^MNS = 


cMNS 

°13 


eMN£ 

~~ °12 


« 90.0° , 






OL\ = 


2(5f 2 NS + 5 2 T S ) = 


2< 2 NS 






02 = 


25 2 M 3 NS = 0° , 







(3.19) 



'23 



9.3° , (3.20) 

where ot\ and ct2 are the Majorana phases as in the PDG parameterization where they are contained 
in a diagonal matrix diag(e iai ' 2 ,e 1Q2 ' 2 , 1). 

We note that with the mixing pattern of our model, i.e. tri-bimaximal mixing produced in the 
neutrino sector, and charged lepton mixing corrections only from #j 2 , the leptonic mixing angles 
and the Dirac CP phase 5mns satisfy the lepton mixing sum rule [241 ] 

#mns _ 0MNS cos ( 5mns ) „ arcs in(l/ v / 3) . (3.21) 

The approximately maximal CP violation, i.e. <5mns ~ 90° > affects that although the charged lepton 
corrections generate 0^3 ° ~ 3.3°, the solar mixing angle remains very close to its tri-bimaximal 
value of arcsin(l/\/3). So far, we have discussed the neutrino mixing parameters at the GUT scale. 
To calculate the low energy values we have to take RG running of the parameters into account. 

3.3 Renormalisation Group Corrections 

For a quasi-degenerate neutrino mass spectrum, RG corrections to the neutrino parameters can 
in principle change the high scale values dramatically. However, as has been discussed for type II 



upgraded seesaw models in [ill] and more generally in [251. l26j|. for small tan/3 and small neutrino 
Yukawa couplings (in our example model they are much smaller than y T ) the corrections to the 
mixing angles and CP phases are under control. Setting the small Majorana phases to zero and 



with <5mns ~ 90° , we can estimate in leading order [25j, [26J 

d#12 ^ Vr „;„/oflMNSw „MNS\2 \ m l + m 2 

I 2 
'sol 



dlnW/^o) ~ 32vr2 Sm( ^ 12 K * 23 } Am 2 (3 ' 22) 



d^ NS 
dhi(/i/Vo) 



(3.23) 



jzjMNS ,,2 /_MNS\2|„- , m |2 i /•„MNS\2|„_ _i_ ~, |2 

dy 23 _ Vt gj n (og MNS ) ( - 12 > \ m2 + m z\ + l g i2 ) \mi + m 3 \ 



dln( M / Mo ) 32tt2 v ^ > Am 2 tm 

13 



where /i is the renormalisation scale. In the case of quasi-degenerate neutrino masses we can further 

use the approximation 7713 ~ m<i ~ m\ = tuq. Integrating these equations approximately with the 

parameters on the right side taken constant and equal to their GUT scale values, and plugging in 

these numbers, we obtain the estimated low energy values of the mixing angles as 

2 
ziMNSi ^ ziMNSi 1 n -1 co m (<i n<\ 

V\2 \mt{m t ) ~ y 12 |A/ GU T+ u - iD fQ 1 eV) 2 ' {o.Aii) 

y 13 \mt(mt) ~ W 13 |M GUT , {d.AO) 

2 

M QUT ± 0-01° m ™° ■ (3.27) 



oMNSi ^ ziMNSi _i_nni° _ 

?23 \m t (m t ) ~ <?23 |M GUT ± U.U1 - - 



In the last equation, the "+" applies for a normal neutrino mass ordering, whereas the "-" applies 
for an inverse mass odering, i.e. the case Am 2 tm < 0. It is important to note that both mass 
orderings can be realised in our model. The strong suppression for the RG running of #^ is 
caused by the particular values of the CP violating phases in our model. For similar reasons the 
running of the CP phases themselves is also suppressed, as can be seen using the analytical results 



m 



25, 26]. In summary, RG corrections are under control in our setup and only cause comparatively 
small corrections to the mixing parameters in the lepton sector. 

In summary, the leptonic mixing parameters in our model are compatible with the experimental 
la ranges at low energy which are: 0^ NS = (34.5±1.0)°, 6^ m = (5.7±|;g)° and 6^ m = (42.3±|jj) , 
taken from [32J ]. as long as mo is smaller than the cosmological bounds suggest, mo < 0.2 eV |2|. 

The values for the leptonic mixing angles and Dirac CP phase <5mns > resulting from our assumed 
vacuum alignment and stated in Eqs. ()3.19p and (I3.20p can be tested accurately by ongoing and 
future precision neutrino oscillation experiments [28(. 

3.4 Predictions for Beta Decay Experiments 

The effective mass relevant for neutrinoless double beta decay is 

m ee = |mic 2 2 c 2 3 e iai + m 2 s\ 2 c\ z i^ + m 3 s 2 3 e 2i5MNS | , (3.28) 

while the kinematic mass accessible in the single beta decay end-point experiment KATRIN is 

2, 2 2 2 i 2 2 2 < 2 2 /o nn\ 

m /3 = m i c nCi3 + m 2 s 12 c r3 + m 3 s 13 . (3.29) 

For quasi-degenerate neutrino mass spectrum (mo = mi — ini — WI3) we obtain that 

m/3 ~ TUo (3.30) 

directly gives information about the absolute neutrino mass scale. 

On the other hand, due to the phases appearing in Eq. (|3.28p there is typically a sizable 
ambiguity in the relation between m ee and ttt-o, as long as the Majorana CP phases are not predicted. 
Allowing, for instance, for arbitrary Majorana phases and considering a quasi-degenerate neutrino 
mass spectrum (mo = mi ~ 777,2 — ^3) an d with small #i3 NS , m ee can still be in the approximate 
interval m ee G [mutest A mutest]- 

This ambiguity is resolved in our model, since the large contribution to the neutrino mass matrix 
proportional to the unit matrix[j (with |mo| 3> \m 3 \, (m^l) results in small Majorana CP phases 



6 In our model this part of the neutrino mass matrix is induced by a standard type I (with right-handed neutrinos) 
and an additional type II seesaw contribution. Nevertheless, the conclusions remain the same as for the pure type II 
seesaw case since we assume that the type II contribution is dominant. 
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Figure 5: The effective mass m ee in our setup relevant for neutrinoless double beta decay as a function of the mass 
m-iightcst of the lightest neutrino, for an inverted neutrino mass ordering (Am^ < 0, upper line) and for a normal 
mass ordering (Am^ > 0, lower line). The bands represent the experimental uncertainties of the mass squared 
differences. The mass bounds from cosmology 2] and from the Heidelberg-Moscow experiment [29(] are displayed as 
grey shaded regions. The red lines show the expected sensitivities of the GERDA experiment in phase I and II [301 ] - 
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and thus we predict: 

m ee « m • (3.31) 

The assumed dominance of mo in our model allows to realise a quasi-degenerate neutrino spectrum 
(with normal or inverse mass ordering) in a natural way, without any tuning of parameters. The 
possible values for m ee as function of the lightest neutrino mass miightest is shown in Fig. [5) 

We would like to remark that for smaller m® one can also naturally extend the model to 
hierarchical or inverted hierarchical neutrino masses without changing the leptonic mixing angles. 
With |mo| ~ |m,2 1 or l m o| ~ \ m i\ or both we then also encounter cases where the Majorana phases 
are close to it. For a quasi-degenerate spectrum the model disfavours these unnatural cases since 
they would correspond to heavily fine-tuned parameters of the model. Similarly, an inverse strongly 
hierarchical spectrum would require unnatural tuning between mo and m^ to make 1713 = \itiq — ml\ 
very small. By contrast, for a typical parameter choice of the model, a normally ordered hierarchical 
spectrum simply corresponds to |mo| <C I^IIj l m 3l an d does not require any tuning at all. It is also 
interesting to note that with the phases in our model there is no possibility to have cancellations 
in Eq. (|3.28[> that could make m ee vanish exactly! 7 ] Neutrinoless double beta decay is thus, also for 
smaller mo, an unavoidable consequence in this class of models. 

4 Summary and Conclusions 

We have proposed a model of quasi-degenerate neutrinos which predicts the neutrinoless double 
beta decay mass observable to be approximately equal to the neutrino mass scale lightest j thereby 
allowing its determination approximately independently of unknown Majorana phases. In general 
such quasi-degenerate neutrino masses may be naturally realised if there exists a non-Abelian family 
symmetry with real triplet representations that enforces an additional contribution to the neutrino 
mass matrix proportional to the unit matrix, hence determining the neutrino mass scale. In our 
model, the additional contribution was generated by a type II seesaw, or, alternatively, by another 
type I seesaw contribution from an additional triplet representation of right-handed neutrinos (or 
neutrino messenger fields). In addition, the standard type I seesaw contribution determines the 
neutrino mixing angles. Although such a mechanism (called a "type II upgrade" in [ll() has been 
known for some time, the model in this paper is the first of its kind to combine this mechanism 
with tri-bimaximal mixing arising from A4 family symmetry together with a SU(5) SUSY GUT. 

The SUSY A4 x SU(5) model considered here has several attractive features. The full renor- 
malisable superpotential for the coupling of matter to flavons, messenger fields and Higgs fields 
is specified, and only discrete auxiliary symmetries are introduced, rather than the more com- 
mon continuous Abelian symmetry that is typically invoked in such models. In the considered 
model, the neutrino mass scale can either originate from the induced vacuum expectation value 
of a 5[/(2)l triplet contained in a fifteen-dimensional representation of SU(5), or from an addi- 
tional contribution to the neutrino mass operator induced by SU(5) singlet messenger fields in 
the triplet representation of A4. Since the type I seesaw contribution involves very hierarchical 
additional neutrino mass contributions, responsible for the tiny mass splittings between the quasi- 
degenerate neutrinos, we use the constrained sequential dominance mechanism which is well suited 
for achieving such strongly hierarchical contributions. In such a model the neutrino flavour sym- 
metry associated with tri-bimaximal mixing is achieved indirectly from the family symmetry, given 
the assumed vacuum alignment. We have included renormalisation group corrections to the mixing 



7 In fact we find numerically that rn ee > 0.007 eV. 
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angles and shown that they are under control for small values of tan/3 in our framework and do 
not significantly modify the leptonic mixing angles. We have therefore considered small values of 
tan/3 throughout our study. 

The model has several interesting phenomenological features which emerge from our numerical 
fit to the quark and lepton masses and quark mixing. In addition to b-r unification that may be 
viable for small tan/3, the model realises the GUT scale relation y^/y s ~ 9/2 proposed in [16], 
which is favourable for small values of tan/3. In the quark sector, we observe that the correct 
<5ckm (corresponding to the right unitarity triangle with a ~ 90°) can be realised for the simple 
ansatz in Eq. (|2.10p . We remark that this simple ansatz leads to an interesting alternative texture 
(different from [22J) that gives rise to a right-angled unitarity triangle with a ~ 90°. For the 
leptonic mixing angles and CP phases we find 6>^ NS » 3.3°, 6»^ NS « 35.1°, 6^ NS » 45.0° and 
<5mns ~ 90° . The leptonic mixing angles satisfy the "lepton mixing sum rule" proposed in 24] with 
only small theoretical errors since the 1-3 mixing in the charged lepton mass matrix is very small 
and the 2-3 mixing vanishes. The Majorana CP phases are small for quasi-degenerate neutrino 
masses via a large additional contribution proportional to the unit matrix, as expected, and the 
model thus predicts the neutrinoless double beta decay mass observable to be approximately equal 
to the neutrino mass scale, or lightest neutrino mass, i.e. m ee ss mutest- 

In conclusion, if neutrinoless double beta decay were observed in the near future then this would 
herald another neutrino revolution in which neutrino masses would be quasi-degenerate. Amongst 
the many possible models of quasi-degenerate neutrinos, the model with an additional contribution 
to neutrino mass matrix proportional to the unit matrix, as considered here, are distinguished by 
their prediction of that the neutrinoless double beta decay mass observable is approximately equal 
to the neutrino mass scale ^lightest- We have proposed the first realistic model of this kind involving 
A4 and SUSY SU(5) GUTs. The A4 family symmetry has the dual effect of enforcing on the one 
hand quasi-degeneracy via the additional unit matrix contribution to the neutrino mass matrix, 
and on the other hand tri-bimaximal mixing via the type I seesaw mechanism and constrained 
sequential dominance. A numerical fit to quark masses and mixing angles reveals that a simple 
ansatz describing quark CP violation with a « 90° also leads to the leptonic phase 5mns ~ 90°. In 
such models the absolute neutrino masses could be directly measurable by experiment quite soon. 
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A The Renormalisable Super potential and Effective Couplings 

With the field content and symmetries specified in Tab. [1] the superpotential contains the following 
renormalisable terms: 

Wh = H5H5H5 + /U15-H15.H15 + [145H45H45 + \15H15H5H5 + \15H15H5H5 , (A.l) 

W A = M Al0 A 10 A 10 + M A5 A 5 A 5 + M Al A\ + M B BB + M Cl (7id + M C2 C 2 C 2 (A.2) 

Wint = KF-F#5^10 + K F F(f)23^5 + KTiTifaAio + KTlTlH^A^ + vaH 15 FF 

r 2 

f 3 



+ A^io^io + Aff 5 C 2J B + /^BT 2 + 4C 2 2 23 + ^C 2 ^ 3 + a t tf 5 T 3 2 



+ RfFH^Ax + « w JVi^Ai + &iQ$ + ZmCiN? . (A.3) 

As discussed in the main text, after GUT symmetry breaking the SU(2)l doublet components 
from H5 and H^ respectively H§ and H^ mix and only the light states acquire the SU(2)l breaking 
vevs which give the fermion masses. We parameterise the Higgs mixing with the mixing angles 7 
respectively 7 

H s\ = fc y -s 7 \ (Hi 

H i5 J " Vs 7 c 7 J \H h/ ( 

#5 \ (Of — ^7 \ / i?Z 



where we have used the common abbreviation c 7 = cos 7 and similar for the sine and the other 
angle 7. The light Higgs doublets are denoted with an index I while the heavy Higgs doublets are 
denoted with an index h. 

The effective couplings a appearing in the effective superpotential in the main text can be 
expressed in term of the fundamental couplings from Eq. ()A.3|) and the messenger masses from 

dA2D 

Oj = KfKTi , a 2 = KFKTi , (A. 5) 



an — Ak^i , a 22 — \Kj<2 + ^ K T K <b~T7 — 77 — ' ° 33 — a < + 



'M B M C2 ' M l M\ 



A w 



Ml 



(A.6) 



a 22 = Xk^ 5 , atj = \K Ti K Tj , for i ^ j , 

r M B M C2 

- Ma w /* v \ 

a Vi = K Ni K F — — , (A.7) 

jyiAi 
2 M\ w M Aw M\ io 

aR t = KNi^TT 21 + ifaiNi-jyjr^- , ClR 12 = KNiKNj^J 21 • (A.8) 

The effective couplings 6 appearing in the Yukawa couplings can be expressed in terms of the 
couplings a and the Higgs mixing angles as 

bi = c^cii , b 2 = s 7 a 2 , (A. 9) 

2 

bij = -f — - , o 22 = -^^2" , 033 = c 7 a 4 + -£-§ , (A. 10) 
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K = ^ , (A.11) 

h - l aRi h - 1 ttRl2 (A 19\ 

B A Possible Vacuum Alignment 

In this section of the appendix, we will discuss a possibility to extend the model in the main part 
by a viable vacuum alignment. This possibility is based on realising the flavon 023 effectively by 
splitting it up into two flavons 02 and 03, which get purely real or purely imaginary vevs by a 
vacuum alignment as described below, i.e.: 

(023> = (02> + (03) where (0 2 ) = I _i I ^23 and (<f> 3 ) = I J e 23 . (B.l) 

W w 

For the effective superpotential in Eqs. (|2.2p - f|2.T[> this would mean to simply replace 023 —> 02 + 03 
and 023-^45 - >• 02-^45 + 03-^5 j with an additional Higgs field H' 5 . The field content and the 
symmetries of this extended version of the model is given in Tab. 01 We note that, at the level of 
precision discussed here, the model predictions are the same as in the main part of the paper. 

We now turn to the discussion of the vacuum alignment sector with the desired purely real or 
imaginary alignment (see also [35(). At the effective theory level the new operators involving, the 
"driving fields" Pi, D{ and Oi, that are generated by the fields and symmetries specified in Tab. [J] 
are given by (dropping 0(1) coupling constants): 

Wz = Pi(Sr-Mn+Difa* fa + OtjMj , (B.2) 




M h 



tA 



W^ = Dl 23 (£0123 + 0123 * 0123) + Pl23 I j^~ + Jr- + ^- - Mf 23 , (B.3) 

\ Cl23 Cl23 Cl23 ) 

Wfc = P 3 (0l - Ml) + D 3 (03 * 03) + 0130103 + O 23 0203 , (B.4) 

Wfea = P23 I jjjg- - M 2 4 3 J + Ol ;23 01023 + Ol23;230123023 ■ (B.5) 

The minimisation of the resulting F-term potential leads to the desired vacuum alignment for 
all flavons. We note that the new flavons £ and 0i are only "auxiliary" in the sense that they do 
not couple to the matter sector but are only relevant for the flavon alignment. The flavon £ leads 
to an additional term contributing to the charm mass, however its effect can be absorbed in the 
messenger masses such that the predictions of the model are unchanged. The four parts Wg, W^ 123 , 
Wfc and W0 23 have the following meaning: 

• The terms in Wg enforce non- vanishing vevs of the three flavons 0i, 02 and 03 that are 
orthogonal to each other and have only one nonzero element, due to the effects of the terms 
with the driving fields Dj. Here the * denotes the symmetric triplet combination of A^ (see 



e.g. [34(]). With real Mj due to the assumed spontaneous CP violation cpj are forced to be real 



which means that the vevs (0j) are either purely real or imaginary. We will choose a vacuum 
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Table 4: Representations and charges of the superfields. The subscript i on the fields Ti and Ni is a family index 
while for Pi and Pi the i matches the corresponding flavon field. The flavon fields fa, $23 can be associated to a 
family via their charges under Z2 x l\. The subscripts on the Higgs fields H, H and extra vector-like matter fields 
A, A denote the transformation properties under SU(5). 
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where 02 is imaginary and 03 is real. The vev of 0i can be either real or imaginary and since 
its phase does not affect the results. 

• The part W0 123 leads to the desired alignment for the flavon vev ((^123) along the (±1, ±1, ±1) 
direction (see also |6|]) in the minima of the potential with (^123) 7^ and (£) 7^ 0. Again, 
due to the 0^ 23 in the term with P123 the vev can be either real or imaginary, and we choose 
the real vacuum in the (1, 1, 1) direction. 

• The terms in W<4 3 and W<£ 23 finally provide the alignment for the flavons 03 and 023- The 
terms with the driving fields 0{ force the vev of 03 to be orthogonal to (0i) and (fa), and 
the vev of 023 to be orthogonal to (0i) and ((^123) , leading to the alignments in the desired 
direction in flavour space. We can again choose the purely real vacuum as explained above. 

We note that the driving fields Pi and Pi have the same quantum numbers and therefore, a 
priori, a linear combination of such fields couples to each of the above terms (ftf/M^, — Mf and 
(J^/Mq. — M23. Without loss of generality we have written the driving fields here already in the 
field basis where only one driving field couples to one of these terms. This can always be achieved 
by a proper field redefinition. 

References 

[1] S. F. King, larXiv:0712. 17501 [physics.pop-ph]; S. F. King, Rept. Prog. Phys. 67 (2004) 



107 



arXiv:hep-ph/0310204|; R. N. Mohapatra et al, Rept. Prog. Phys. 70 (2007) 1757 



arXiv:hep-ph/0510213|; R. N. Mohapatra and A. Y. Smirnov, Ann. Rev. Nucl. Part. Sci. 



56 (2006) 569 |arXiv:hep-ph/0603118|; C. H. Albright, arXiv:0905.0146 [hep-ph]; G. Altarelli 
and F. Feruglio, [arXiv:1002.02H] [hep-ph]. 

[2] N. Jarosik et al. larXiv: 1001 .47441 [astro-ph.CO]. 

[3] F. Feruglio, A. Strumia and F. Vissani, Nucl. Phys. B 637 (2002) 345 [Addendum-ibid. B 659 
(2003) 359] [arXiv:hep-ph/0201291] . 

[4] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530 (2002) 167 
larXiv:hep-ph/0202074] ; P. F. Harrison and W. G. Scott, Phys. Lett. B 557 (2003) 76 
|arXiv:hep-ph /0302025| . 

[5] S. F. King and C. Luhn. larXiv:0912. 13441 [hep-ph], 

[6] G. Altarelli and F. Feruglio, Nucl. Phys. B 741 (2006) 215 [arXiv:hep-ph/0512103] . 

[7] G. Altarelli, F. Feruglio and C. Hagedorn, JHEP 0803 (2008) 052 [arXiv:0802.0"nM1 [hep-ph]]: 
T. J. Burrows and S. F. King. larXiv:0909.1433l [hep-ph], 

[8] P. Ciafaloni, M. Picariello, E. Torrente-Lujan and A. Urbano, Phys. Rev. D 79 (2009) 116010 
[arXiv:0901.2236 [hep-ph]]; P. Ciafaloni, M. Picariello, A. Urbano and E. Torrente-Lujan, Phys. 
Rev. D 81 (2010) 016004 [arXiv:0909.2"553l [hep-ph]]. 

[9] E. Ma and G. Rajasekaran, Phys. Rev. D 64 (2001) 113012 [arXiv:hep-ph/0106291]; E. Ma 



arXiv:hep-ph /0701016; E. Ma, Mod. Phys. Lett. A 22 (2007) 101 |arXiv:hep-p h/0610342 1 ; 



21 



E. Ma, Mod. Phys. Lett. A 21 (2006) 2931 |arXiv:hep- ph/0607190 j ; E. Ma, H. Sawanaka 
and M. Tanimoto, Phys. Lett. B 641 (2006) 301 [arXiv:h~p h/0606103|; E. Ma, Phys. Rev. 
D 73 (2006) 057304; B. Adhikary, B. Brahmachari, A. Ghosal, E. Ma and M. K. Parida, 
Phys. Lett. B 63 8 (2006) 345 [arXiv:hep-ph7 0603059 1 ; E. Ma, Mod. Phys. Lett. A 20 (2005 ) 



2601 arXiv:hep-ph/0508099|; E. Ma, Phys. Rev. D 72 (2005) 037301 [arXiv:hep-ph/0505209] 



S. L. Chen, M. Frigerio and E. Ma, Nucl. Phys. B 724 (2005) 423 |arXiv:hep-ph/0504181|; 
W. Rodejohann and Z. z. Xing, Phys. Lett. B 601 (2004) 176 |arXiv :hep-ph/0 408195]. E. Ma, 
Phys. Rev. D 70 (2004) 031901 |arXi v:hep-ph/0404199|; F. Feru glio, C. Hagedorn, Y. Lin and 
L. Merlo, Nucl. Phys. B 775 (20 07) 120 [arXiv:hep-ph/0702194| ; X. G. He, Nucl. Phys. Proc. 
Suppl. 168 (2007) 350 |arXiv :hep-ph/0612080| ; S. F. King and M. Malinsky, Phys. Lett. B 645 
(2007) 351 |arXi v:hep-ph/0610250] ; I. K. Cooper, S. F. King and C. Luhn, arXiv: 1004.3 243 
[hep-ph] . 

[10] R. Barbieri, L. J. Hall, G. L. Kane and G. G. Ross, |arXiv:hep-ph7 9901228 ; I. de Medeiros 
Varzielas and G. G. Ross, Nucl. Phys. B 733 (2006) 31 |arXiv:hep-ph/0507176|; S. F. King 



and M. Malinsky, JHEP 0611 (2006) 071 |arXiv:hep-ph/0608021 ; I. de Medeiros Varzielas 



S. F. King and G. G. Ross, Phys. Lett. B 644 (2007) 153 |arXiv:hep-ph/0512313|; 
I. de Medeiros Varziel as, S. F. King and G. G. Ross, Phys. L ett. B 648 (2007) 201 
|arXiv:hep-ph /0607045 j ; E. Ma, Mod. Phys. Lett. A 21 (2006) 1917 |arXiv:hep-ph/0607056|; 
E. Ma, larXl v:0709.0507 [hep-ph]; C. Luhn, S. Nasri and P. Ramond, Phys. Lett. B 652 
(2007) 27 [arXiv:0706.2341 [hep-ph]]; F. Plentinger and G. Seidl, Phys. Rev. D 78 (2008) 
045004 [arXiv:0803.2889 [hep-ph]]; C. Csaki, C. Delaunay, C. Grojean and Y. Grossman, 
arXiv:0806.0356 [hep-ph]; M.-C. Chen and K. T. Mahanthappa, Phys. Lett. B 652 (2007) 34 
|arXiv:0705.0"714l [hep-ph]]; M.-C. Chen and K. T. Mahanthappa, larXiv:0710.2TT8l [hep-ph]; 
M.-C. Chen and K. T. Mahanthappa, arXiv:0812.4981 [hep-ph]; R. N. Mohapatra, S. Nasri 
and H. B. Yu, Phys. Lett. B 639 (2006) 318 |arXiv:hep-ph/0605020J; R. N. Mohapatra and 
H. B. Yu, Phys. Lett. B 644 (2007) 346 [arXiv:hep-p h/0610023|; A. Aranda. larXiv707 07.3661 
[hep-ph]; A. H. Chan, H. Fritzsch and Z. z. Xing. larXiv:0704.3153l [hep-ph]: Z. z. Xing, Phys. 
Lett. B 618 (2005) 141 [arXiv:hep- ph/0503200 1 ; Z. z. Xing, H. Zhang and S. Zhou, Phys. 
Lett. B 641 (2006) 189 |arXiv:hep-ph/0607091]; S. K. Kang, Z. z. Xing and S. Zhou, Phys. 
Rev. D 73 (2006) 13001 [arXiv:hep-ph/0511157] ; S. Luo and Z. z. Xing, Phys. Lett. B 632 
(2006) 341 | arX iv:hep-p h/0509065 1 ; M. Hirsch, E. Ma, J. C. Romao, J. W. F. Valle and 
A. Villanova del Moral, Phys. Rev. D 75 (2007) 053006 [arXiv:hep-ph /0606082 1 ; N. N. Singh, 
M. Rajkhowa and A. Bora h, |arXiv:hep- ph/0603189; X. G. He and A. Zee, Phys. Lett. B 
645 (2007) 427 [arXiv:hep-p h/0607163 ] ; N. Haba, A. Watanabe and K. Yoshioka, Phys. 
Re v. Lett. 97 (2006) 0416 01 [arXiv:hep-ph/0603116] ; Z. z. Xing, Phys. Lett. B 533 (2002) 
85 |arXiv:hep-ph /0204049 1 ; Y. Lin, Nucl. Phys. B 813 (2009) 91 [arXiv:0804.2867 [hep-ph]]; 
L. Yin, arXiv:0903.0831 [hep-ph]; S. F. King and C. Luhn, arXiv:0905.1686 [hep-ph]; C. Hage- 
dorn, S. F. King and C. Luhn, arXiv: 1003.4249 [hep-ph]. 



[11] S. Antusch and S. F. King, Nucl. Phys. B 705 (2005) 239 |arXiv:hep-ph/0402121| 



[12] S. F. King, Phys. Lett. B 439 (1998) 350 |arXiv:hep-ph/9806440]; S. F. King, Nucl. 
Phys. B 562 (1999) 57 |arXiv:hep-ph/9904210|; S. F. King, Nucl. Phys. B 576 (2000) 85 
|arXiv:hep-ph/9912492|; S. F. King, JHEP 0209 (2002) 011 |arXiv:hep-ph/0204360|; S. An- 
tusch and S. F. King, New J. Phys. 6 (2004) 110 |arXi v:hep-ph/0405272| ; S. Antusch, S. Boud- 
jemaa and S. F. King. larXiv:T003 .5498 [hep-ph]. 



22 



[13] S. F. King, JHEP 0508, 105 (2005) |arXiv:hep-ph/0506297|. 



[14] S. F. King and C. Luhn, JHEP 0910 (2009) 093 [arXiy:0908.1897 [hep-ph]]. 

[15] I. d. M. Varzielas, G. G. Ross and M. Serna, Phys. Rev. D 80 (2009) 073002 [arXiv:081 1.22261 
[hep-ph]]. 

[16] S. Antusch and M. Spinrath, Phys. Rev. D 79 (2009) 095004 [arXiv: 0902 .46441 [hep-ph]]. 



[17] G. Ross and M. Serna, Phys. Lett. B 664 (2008) 97 [arXiv:0704.1248 [hep-ph]]. 



[18] R. Howl and S. F. King, arXiv:0908.2067 [hep-ph]. 

[19] B. Bajc, P. Fileviez Perez and G. Senjanovic, |arXiv:hep-ph/0210374J 



[20] S. Antusch and M. Spinrath, Phys. Rev. D 78 (2008) 075020 |arXiv:0804.0717l [hep-ph]]. 
[21] A. E. Nelson, Phys. Lett. B 136 (1984) 387; S. M. Barr, Phys. Rev. D 30 (1984) 1805. 
[22] S. Antusch, S. F. King, M. Malinsky and M. Spinrath. larXiv:0910.5127l [hep-ph]. 
[23] Z. z. Xing, H. Zhang and S. Zhou, Phys. Rev. D 77 (2008) 113016 [arXiv:0712.14l9l [hep-ph]] . 



[24] S. F. King, JHEP 0508 (2005) 105 |arXiv:hep-ph/0506297|; I. Masina, Phys. Lett. B 633 
(2006) 134 |arXiv:hep-p h/0508031 1; S. Antusch and S. F. King, Phys. Lett. B 631 (2005) 42 
|arXiv:hep-ph/0508044|; S. Antusch, P. Huber, S. F. King and T. Schwetz, JHEP 0704 (2007) 



060 arXiv:hep-ph/0702286|. 



[25] S. Antusch, J. Kersten, M. Lindner and M. Ratz, Nucl. Phys. B 674 (2003) 401 
|arXiv:hep-ph70 305273| . 

[26] S. Antusch, J. Kersten, M. Lindner, M. Ratz and M. A. Schmidt, JHEP 0503 (2005) 024 
[arXiv:hep-ph/0501272] . 

[27] S. Antusch, S. F. King and M. Malinsky, Nucl. Phys. B 820 (2009) 32 [ arXiv:0810.3863l [hep- 
ph]]; S. Antusch, S. F. King and M. Malinsky, JHEP 0805 (2008) 066 [ arXiv:0712.3759l [hep- 
phjj; S. Antusch, S. F. King and M. Malinsky, Phys. Lett. B 671 (2009) 263 [arXiv:0711.4727l 
[hep-ph]]; S. Boudjemaa and S. F. King, arXiv:0808.2782l [hep-ph]; S. F. King, Phys. Lett. B 
659 (2008) 244 [ arXiv:0710.0530l [hep-ph]]. 

[28] A. Bandyopadhyay et al. [ISS Physics Working Group], Rept. Prog. Phys. 72 (2009) 106201 
|arXiv:0710.4947l [hep-ph]]. 

[29] H. V. Klapdor-Kleingrothaus et al, Eur. Phys. J. A 12 (2001) 147 | ar Xiv:hep-ph/0103062| . 

[30] A. A. Smolnikov [GERDA Collaboration], [arXjyiOJilJL4194] [nucl-ex]. 



[31] A. Osipowicz et al. [KATRIN Collaboration], |arXiv:hep-ex/ 0109033 , 

[32] M. C. Gonzalez-Garcia, M. Maltoni and J. Salvado. larXiv: 1001 .45241 [hep-ph]. 

[33] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 (2008) 1. 



[34] S. F. King and M. Malinsky, Phys. Lett. B 645 (2007) 351 |arXiv:hep-ph/0610250 



[35] S. Antusch, S. F. King, C. Luhn and M. Spinrath. lar"Xiv:1103.5930l [hep-ph]. 

23 



